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Abstract 

The Lattice results for the pion and kaon decay constants are analysed within the 
Resonance Chiral Theory framework in the large Nq limit. The approximately linear 
behaviour of the observable at large light-quark mass is explained through the interac- 
tion with the lightest multiplet of scalar resonances. The analysis of the Lattice results 
allows to obtain the resonance mass Ms = 1049 ± 25 MeV and the Chiral Perturbation 
Theory parameters at leading order in 1/Nc- 



1 Introduction 



Quantum Chromodynamics (QCD) has been proved to be the proper theory to describe 
the strong interactions. However in the low energy region the theory in terms of quarks 
and gluons becomes highly non perturbative. These degrees of freedom get confined within 
complex hadronic states. Below the first resonance multiplet (E <C M p ), the only particles 
in the spectrum are the light octet of pseudo-scalars, the pseudo-Nambu-Goldstone bosons 
(pNGB) from the spontaneous chiral symmetry breaking. At low momenta and small pNGB 
masses one may describe their interactions through a chiral invariant effective field theory, 
Chiral Perturbation Theory (xPT) |TJ . It establishes an expansion on powers of the external 
momenta and masses over a characteristic chiral scale A x ~ 4ttF ~ 1.2 GeV, being F ~ 
F n = 92.4 MeV the physical pion decay constant. 

The chiral expansion breaks down when either the momenta or the pNGB masses be- 
come large; when they approach to the p(770) mass the chiral expansion fails. The mesonic 
resonances can then be produced and their effects cannot be neglected any longer. 

Alternatively, it is possible to employ the l/Nc expansion to describe the matrix elements 
-being Nc the number of colours- Jjj. Resonance Chiral Theory (RxT) [3| incorporates the 
interactions between the resonances and the Goldstones at leading order in l/Nc and also 
implements the chiral symmetry of the interaction. Likewise, ;\RT is fully recovered when 
going to low energies. 

At leading order in l/N c (LO) the observables are given by the tree- level amplitudes, 
being the mesonic loops suppressed by l/Nc |2]. In that situation R%T is able to reproduce 
the short distance behaviour required by QCD for the pion form factors, two-point Green 
functions and forward scattering amplitudes [HE]. The l/Nc counting can also be carried 
to the next order (NLO) in a systematic way and the quantum loop corrections might be 
calculated 

The present study is focused on how the variation of the quark masses affects the pion and 
kaon decay constants, F n and Fk, under the R^T framework and the the l/Nc expansion. 
Lattice calculations have provided information about QCD results for unphysical values of 
the u/d light quarks 7 . The simulations are forced sometimes to work with masses of nearly 
the size of the physical strange quark mass or higher. Thus, the usual %PT extrapolations 
break down and generate large unphysical chiral logarithms, yielding an large bending in 
the extrapolation [HEH^]. Nonetheless, more than purely numerical values, this work aims 
providing a possible way to analyse the Lattice simulations at large quark masses, explain- 
ing why the usual linear extrapolations work so well and what are the underlying physical 
foundations. 

Through the inclusion of the first resonance multiplets, with masses Mr ~ 1 GeV, one 
expects to reproduce the physics for the Goldstones up to that range of momenta and masses. 
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Moreover, at LO in 1/Nc only the scalar resonances contribute to F n and Fk- Since the first 
multiplet of pseudo-scalar resonances is at much higher masses, it will not be considered in 
the calculation and its mixing with the pNGB will be neglected. 



2 RxT lagrangian at LO in 1/N ( 



c 

The leading order resonance lagrangian was developed in Ref. [Sj and includes one multi- 
plets of vector, axial- vector, scalar and pseudo-scalar resonances. It is composed by a pair of 
terms including only Goldstones -the 0{p 2 ) x?T Lagrangian pQ , 

F 2 

A = T («X + X+) , (1) 

and pieces including as well resonance fields. The brackets (...) mean trace of the flavour 
matrices. The chiral tensors including the Goldstone fields can be obtained in Refs. j3 IH El 
and they contain the external fields v^, a M , s and p, and the non-linear realization of the chiral 
symmetry given by the tensor u = exp (i$/V2F), with the pNGB fields $ = Ea<t>aK/ y/2. 

In the large Nq limit the qq resonances form U(3) multiplets. The fields of a multiplet 
can be put together in a 3 x 3 matrix which transforms linearly under the chiral symmetry, 

containing one chiral singlet field i?o and the remaining octet fields R a , with a — 1, ...8. All 
the fields FL a of the multiplet would be degenerate in the large N c and massless quark limit. 
The kinetic terms are then constructed with these matrices: 



(3) 



Cf a (R = V,A) = -\{V x R Xf N v R^ -IMIR^R^) , 

Cf»(R = S,P) = \ ( V^RV^R - M 2 R R 2 } , 
and the interaction lagrangian, linear in the resonance fields jH], 

C ' 2A - &<^>- (4) 

C-2S = C d (Su^) + C m (SX+) , 

C-2P = id m (PX-) > 

where the vector and axial-vector resonances are described by antisymmetric tensors. 



2 



3 Scalar tadpole and field redefinition 



Analysing the LO lagrangian one observes the presence of a term linear in the scalar fields, 
i.e. a scalar tadpole. It is given by the term of C 2 s in Eq. (jU with the coupling c m , which 
also provides the vertex for the scalar resonance production from a scalar quark current. 

Chiral symmetry requires the quark masses to enter in the effective lagrangian only 
through the tensor \ = 2Bq{s(x) + ip(x)}, which appears in the chiral covariant com- 
binations x± = u^xu) ± uy}u. In order to recover physical QCD, the external fields are 
evaluated at the end of the calculation as x = 2B A4, being M. =diag(m u , m d , m s ) the 
diagonal matrix with the light quark masses. 

The pieces of the lagrangian containing the scalar fields are 

£f«(S) + C 2S = \(d»Sd,S) - \m 2 {S 2 > + 4B c m (SM ) + 0(S 2 ® 2 , S$ 2 ) . 

(5) 

The scalar field has therefore a non-zero vacuum expectation value (v.e.v.). In order to define 
the quantum field theory around the minimum one needs to perform in the scalar field the 
shift: 

+ (6) 

where the shifted scalar nonet fields S has a zero v.e.v. Nonetheless this shift makes S not to 
be chiral covariant any longer. For the present work this detail is not relevant although other 
alternative shifts (like S = S + c m x+) would restore the explicit covariance. The important 
detail is that the shift is not equal for all the scalar fields of the multiplet, but proportional 
to the quark mass matrix M. and different for each resonance. 

The part of the lagrangian containing the vectors, axials and pseudo-scalar resonances, 
Cr[V, A, P], remains unchanged under the shift but the remaining 0(p 2 ) chiral term C 2x and 
the scalar pieces £f in + C 2 s become 

Cf« + C 2S + C 2x = Cf-' + C' 2S + C' 2x - ^(-M 2 ) , (7) 

yielding a constant term proportional to (Ai 2 ), a kinetic term structure for S, 

4 m ' = ~ (V- (s + ^fM^j V, (s + ^M^j ) - \ ( MIS 2 ) , (8) 

an interaction lagrangian without tadpoles, 

C' 2 - s = c d ( Su^ ) + c m ( S ( X + ~ 4B M) ) (9) 
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and a modified 0(p 2 ) %PT term, 



r - — 

~ 4 

provided by the matrix definitions 



( cr 1 k«" + x + ] ) 



/^-i | i i 4c d c m 4£ n A'i 



-i 



1 , 4cj 4gpAj 

1 + iTS - 



F z M; 
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(10) 



In order to convert the pNGB kinetic term to the canonical form, one needs to perform 
a re-scaling on the pNGB fields u = exp (i$/y/2F) = exp (iC^^ /y/2F\ . 



At LO in l/iVc the resonance couplings are fixed by the QCD short distance constraints 5, 
ITU] : Cd = c m = F/2. Hence, from Eqs. (JTUJ) and (JTJJ) at LO in 1/Nc one gets for the pion 
and kaon fields the re-scalings and masses 



ml 



2B m . 



2m 
1 4_ V 



m 



K 



Mi 



C 



K 



B (m + m s ) 



1+ Ml" 



(12) 



3.1 and F# at Leading Order 



Since the scalar tadpole has been removed, at LO there is only one diagram contributing 
to the pion decay constant: the tree-level production of the pNGB from the axial current, 



0\d lfll5 u\P + (p)) 



V2FC 



and therefore the pion and kaon decay constants are 

2m! 



F 1 + 



Mi 



K 



p 2 V>i 



F 1 + 



2m 2 K y 



Mi 



(13) 



(14) 



When mp <C M| the decay constants may be expanded in powers of Trip, recovering the 



tree-level xPT result F P = F 1 + ^mp + ^{2m 2 K + 



ml 



, with L 5 = F 2 /AM 2 
1. Only near the zero 



0{m%) 

and L 4 = [5]. This explains why the linear extrapolations work so we 
quark mass the LO result in 1/Nc gains sizable non-analytic contributions from the one- loop 
logarithms (NLO in 1/N C ): F n = 



1 + 



F 2 



ml 



16vr 2 F 2 



In =* + ... 
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Figure 1: A) Lattice results for F n and Fk and extrapolations taken from MILC collabora- 
tion 0. B) Comparison of the Lattice results with the R%T extrapolations for the values from 
the fit, F = 93.7 MeV, M s = 1020 MeV and the splitting e s m = -0.02. The kaon mass varies as 
m 2 K = (m K hys ) 2 + [m 2 — (m^ hys ) 2 ]/2 due to its light quark content. They are shown together with 
the data for = rn s /2.3 (squares), = m s /4.5 (triangles) and the linear extrapolation to 

the value = (error bars). 



The results for F n and Fk from the Lattice simulations (MILC Collaboration [7J are 
shown in Fig. (l.A). The simulation handles two kinds of quark masses: The sea-quark 
masses of the fermions within closed loops; and the valence-quark masses of those which are 
not from the sea. In this simulation the strange quark valence-mass m v s al and the strange 
quark sea-mass m s s ea take the same value m s = m v s al = m s s ea = 1.14 mf hys , being mf hys 
the physical mass of the strange quark. It may be related with the physical kaon mass 
through Eq. (|12j). The u/d quark valence- masses (isospin limit is assumed) are varied 
continuously between nearly zero and the physical mass of the strange quark. Finally, the 
simulation is run for two values of the sea-masses of the quarks u/d: = m s /2.3 and 

m t?d = n^s/^-5 (squares and triangles respectively in Fig. (l.B)). 

The modifications due to the sea quark mass are much smaller than those from the 
valence masses, as it is expected in the large Nq limit, since the closed quark loops would be 
suppressed by 1/Nq. Thus, in this work I have generated the matrix elements for equal values 
of the sea and valence masses, i.e. for = m^, through a simple linear extrapolation 
from the samples = m s /2.3 and = m s /4:.5 (error bars in Fig. (l.B)). 

The theoretical expressions derived from R%T for F n and Fk in Eqs. (114)1 were fitted 
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to these extrapolated points. The error in my input data in the fit was a 1 % error, a 
typical discretization error, but it did not account the uncertainties in the extrapolation 
to the value = m™\. The variation of the decay constant when changing m s ^ could 
be considered as a rough estimate of this uncertainty. In addition, the errors due to NLO 
contributions in l/Nc have not been considered. The fit yields the values F = 94.1±0.9 MeV 
and Mg = 1049 ± 25 MeV, with x 2 /dof= 11.0/13. This yields for the physical pion and kaon 
decay constants F n = 95.8 ±0.9 MeV and Fx = 113 ±1.4 MeV, in acceptable agreement with 
the experimental values F w + = 92.4 ±0.07 ±0.3 MeV and F K + = 113.0 ± 1.0 ± 0.3 MeV [TT]. 
One may also estimate the %PT coupling at LO in 1/N C , L 5 = F 2 /AM 2 S = (2.01±0.10) ■ 10~ 3 . 

Nevertheless, here the chiral logarithms are lacking. In they produce an important 
non-analytic effect and a large bending in the F n curve, since its slope becomes large at small 
pion mass due to J^f- ~ In — — . The effect of the logs on Fj^ is much more reduced since 
the value of the kaon mass does not become small when m U) d — > 0. That is the reason for 
the better agreement of the Fk result. Eventually the NLO calculation in l/Nc (one loop) 
would introduce this extra non-analytic curvature and the usual one loop result for F n in 
%PT would be recovered. 

The mass splitting between the masses of the two 1 = scalar resonances can also be 
studied since they do not contribute equally to F w and F K . The discussion is developed 
in the appendix. At LO in l/Nc there would be two mass eigenstates with square masses 
M| n = Mg — 8e^B m Ut d and Mg s = Mj — 8e^-Bo^ s , with quark contents (-t^uu + -j^dd) and 
ss respectively. The pion and kaon decay constants depend now on the splitting parameter e^. 
However the fit to the former data is not sensitive to this coupling, yielding: F = 93.7 ± 1.5 
MeV, M s = 1020±80 MeV and e s m = -0.02±0.05, with similar xV dof = 10.9/12. It provides 
an estimate of the scalar masses at large Nc'- one gets the values Mg n = 1020 ± 80 MeV and 
M Ss = 1040 ± 90 MeV, highly correlated, and their splitting, equal to M Ss - M Sn = 20 ± 40 
MeV. In addition, large OZI-rule violations could occur at NLO of the same size as the quark 
mass corrections. In the worst case, one would expect them to be of the order of l/Nc — 33% 
times the value of the scalar mass at LO in l/Nc- 

The R%T extrapolation is shown in Fig. (l.B) (solid line). Since this effective field theory 
is for equal valence and sea masses, the theoretical result is compared with an emulation of 
the Lattice data for = m^ al d , obtained by linear extrapolation as it was explained before. 



4 Conclusions 

R%T has been shown to be an interesting tool to analyse the Lattice data, which usually are 
generated for non-physical values of the light quark masses. The present study hints that the 
light mesonic resonances may play an important role in the large quark mass extrapolations. 
This work explores the pion and kaon decay constants, providing successful results. Its 
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importance and aim is not just the determination of the decay constants but to present an 
alternative idea to interpret the Lattice simulations for large unphysical values of the masses, 
giving a clear explanation of why the usual linear extrapolations yield such a good result. 
Likewise, this provides solid theoretical foundations based on the underlying QCD for these 
techniques. Thus, the 1/Nc expansion might be as well a suitable framework to describe the 
heavy quark matrix elements Bb---) at large values of the u/d quark masses, where a 
similar linear behaviour has been also observed. 

The fact that at low energies R%T recovers ensures that we are introducing the 

proper low mass behaviour ISHUEI- Former works jHUHj noticed the necessity of a separation 
scale A cut where the xPT loops become irrelevant. The resonance masses provide a "natural" 
scale where the chiral extrapolations fails and where the dynamics of the observable changes 
drastically. 

The fits to the simulations were done for an emulation of the Lattice data, obtained by 
extrapolating m s ^ to the value = m^, that varied in a wide and continuous range 
between zero and the strange quark mass. For a more proper analysis one would need 
a simulation with equal sea and valence masses. However, the main dependence comes 
from the valence-quarks and the sea-quark effects are small, since the closed quark loops 
are suppressed by 1/Nc- Therefore, the present calculation can be considered an adequate 
estimate of the hadronic parameters. From the former fit the values F = 94.1 ± 0.9 MeV 
and Ms = 1049 ± 25 MeV were obtained, together with the xPT coupling estimate at LO in 
1/Nc, L 5 = (2.01 ± 0.10) ■ 10 -3 . The scalar mass splitting showed large uncertainties. 
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Appendix: Scalar mass splitting 

The resonance multiplet has been considered at first as degenerate in mass. Nonetheless 
the resonance masses can gain contributions due to the quark masses. In the large Nc limit, 
Chiral invariance requires that at order 0(m q ) the mass splitting comes only through a chiral 
invariant term [T2j : 

AC mq = e s m ( X+ S 2 ), (15) 
being an 0(N C ) dimensionless constant, independent of the quark masses. 
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The shift now in the scalar field to re-absorb the tadpole is slightly different: 

-l 



S = S + AB c m M M' z - %e s m B Q M 



(16) 



The mass eigenvalues for the 1 = scalars are not M| anymore but the two values 
M Sn = M| — 8e^-B m Uirf , for the state C^uu + ^dd), and M Sa = Mj — 8e^-B ms, for ss. 
Nonetheless, the physical scalar states will separate away from this ideal mixing and these 
masses will gain also contributions due to NLO effects in 1/Nc- 

At LO in 1/Nc the pNGB masses still remain as given in Eq. (|12j). However the re-scaling 
factors for pions and kaons change accordingly: 



1 + 



2m 



Ml - Ae s m ml 



C 



K 



1 + 



2m 



K 



Mi 



4e s m m 2 K 



(17) 



which therefore modify the pion and kaon decay constants. 
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